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Abstract 

A dispersion integral is derived that allows one to relate directly (spin dependent) 
AN invariant mass spectra, measured in a large-momentum transfer reaction such 
as pp — > K~^pA or 'jd — > K~^nA, to the scattering length for elastic AA^ scattering. 
The involved systematic uncertainties are estimated to be smaller than 0.3 fm. This 
estimate is confirmed by comparing results of the proposed formalism with those of 
microscopic model calculations. We also show, for the specific reaction pp — > K~^Ap, 
how polarization observables can be used to separate the two spin states of the AA^ 
system. 



1 Introduction 



Nucleons and hyperons form a flavor-S'[/(3) octet. The underlying SU{3) 
symmetry is clearly broken, as is evidenced by the mass splittings within the 
members of the octet. This symmetry breaking can be well accounted for with 
relations such as the Gell-Mann-Okubo formula. The interesting question is, 
however, whether there is also a dynamical breaking of the SU{3) symmetry 
- besides these obvious "kinematic" effects. 

The nucleon-nucleon {NN) interaction has been studied in great detail for 
many decades and we have a good understanding of this system at energies 
even beyond the pion production threshold. On the other hand, very little is 
known about the dynamics involving the other members of the octet. Because 
of that it has become common practice in studies of hyperon physics to assume 
a priori SU{3) flavor symmetry. Specifically, in meson exchange models of the 
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hyperon-nucleon (YN) interaction such symmetry requirements provide rela- 
tions between coupling constants of mesons of a given multiplct to the baryon 
current, which greatly reduce the number of free model parameters. Then cou- 
pling constants at the strange vertices are connected to nucleon-nucleon-meson 
coupling constants, which in turn are constrained by the wealth of empirical 
information on NN scattering [1,2,3,4,5,6]. The scarce and not very accurate 
data set available so far for clastic YN scattering [7,9,10] seems to be in- 
deed consistent with the assumption of SU{3) symmetry. Unfortunately, the 
short lifetime of the hyperons hinders high precision scattering experiments 
at low energies and, therefore, has so far precluded a more thorough test of 
the validity of SU (3) flavor symmetry. 

The poor status of our information on the YN interaction is most obviously 
reflected in the present knowledge of the AN scattering lengths. Attempts 
in the 1960's to pin down the low energy parameters for the S-wecves led to 
results that were afflicted by rather large uncertainties [7,8]. In Ref. [8] the 
following values are given for the singlet scattering length and the triplet 
scattering length at 

as = -1.8 I fm and at = -1.6 I^Jg fm, (1) 

where the errors are strongly correlated. The situation of the corresponding 
effective ranges is even worse: for both spin states values between and 16 
fm are allowed by the data. Later, the application of microscopic models for 
the extrapolation of the data to the threshold, was hardly more successful. 
For example, in Ref. [5] one can find six different models that equally well 
describe the available data but whose (S'-wave) scattering lengths range from 
-0.7 to -2.6 fm in the singlet channel and from -1.7 to -2.15 fm in the triplet 
channel. 

A natural alternative to scattering experiments are studies of production re- 
actions. In Ref. [11] it was suggested to use the reaction K~d nA'-f, where 
the initial state is in an atomic bound state, to determine the AA^ scattering 
lengths. From the experimental side so far, a feasibility study was performed 
which demonstrated that a separation of background and signal is possible 
[12]. The reaction K~d — > nA'j was studied theoretically in more detail in 
Refs. [13,14,15]. The main results especially of the last work are that it is 
indeed possible to use the radiative capture to extract the AA^ scattering 
lenghts and that polarization observables could be used to disentangle the 
different spin states. In that paper it was also shown, however, that to some 
extend the extraction is sensitive to the short range behavior of the YN in- 
teraction. The reaction K~d Ti~pA was analyzed in Ref. [16], leading to 
a scattering length of —2 ± 0.5 fm via fitting the invariant mass distribution 
to an effective range expansion — the author argued that this value is to be 
interpreted as the spin triplet scattering length. It is difficult to estimate the 
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theoretical uncertainty and the error given is that of the experiment only. 

In the present paper we argue that large-momentum transfer reactions such 
as pp K+pA [17,18,19] or 'yd K+nA [20,21,22,23,24] might be the best 
candidates for extracting informations about the AN scattering lengths. In 
reactions with large momentum transfer the production process is necessarily 
of short-ranged nature. As a consequence the results are basically insensitive to 
details of the production mechanism and therefore a reliable error estimation 
is possible. In Ref. [19] the reaction pp — > K^pA at low excess energies was 
already used to determine the low energy parameters of the AN interaction. 
The authors extracted an average value of —2 ± 0.2 fm for the AA^ scattering 
length in an analysis that utilizes the effective range expansion. But also this 
work has some drawbacks. First, again the given error is statistical only. More 
serious, however, is the use of the effective range expansion. As we will show 
below this is only appropriate for systems in which the scattering length is 
significantly larger than the effective range. In addition, using this procedure 
one encounters strong correlations between the effective range parameters a 
and r that can only be disentangled by including other data, e.g. AA^ elastic 
cross sections, into the analysis [19]. 

In this manuscript we propose a method that allows extraction of the YN scat- 
tering lengths from the production data directly. In particular, we derive an 
integral representation for the AA^ scattering lengths in terms of a differential 
cross section of reactions with large momentum transfer such as pp — > K'^pA 
or 7^ — > K'^nA. It reads 



where as denotes the spin cross section for the production of a A-nucleon 
pair with invariant mass m' ^ — corresponding to a relative momentum p' — and 
total spin S. In addition t = {pi —pK+Yi with pi being the beam momentum, 
mQ = (mA -f-mAr)^, where mA {itln) denotes the mass of the Lambda hyperon 
(nucleon), and nimax is some suitably chosen cutoff in the mass integration. We 
will argue below that it is sufficient to include relative energies of the final AA" 
system of at most 40 MeV in the range of integration to get accurate results. 
P denotes that the principal value of the integral is to be used and the limit 
has to be taken from above. This formula should enable determination of the 
scattering lengths to a theoretical uncertainty of at most 0.3 fm. Note that, 
due to recent progress in accelerator technology, the differential cross section 
required can be measured to high accuracy even for small AA" energies. 
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The method we propose is apphcable to all large momentum transfer reactions, 
as long as the final AA^ system is dominated by a single partial wave. Evidently, 
the number of contributing partial waves is already strongly reduced by the 
kinematical requirements specified above. To be concrete, we do not expect 
that P or higher partial waves are of relevance for KN energies below the 
mentioned 40 MeV. However, in an unpolarized measurement both the spin 
triplet (^5*1) as well as the spin singlet {}Sq) final state can contribute with a 
priory unknown relative weight. Fortunately polarization observables allow to 
disentangle model independently the different spin states. In appendix B we 
demonstrate this for the specific reaction pp — > pAK'^. 



2 Formalism 



2.1 The Enhancement Factor 



A standard method to calculate the effect of a particular final state interaction 
is that of dispersion relations [25]. The production amplitude A depends on 
the total energy squared s — (pi +^2)^, the invariant mass squared of the 
outgoing Ap system = (Pn + PaY and the momentum transfer t defined 
above. The dispersion relation for A at fixed s and t and a specific partial 
wave S takes the form 



As[s,t,m)^- —dm +-I 7^ —dm , (3) 

TT J -00 m' ^ — TT Jml m' ^ — m^ 

where m^ = m ^ is the lefthand singularity closest to the physical region and 



Ds{s, t, m^) = ^(^5(5, t, m^ + iO) - As{s, t, m^ - iO)). (4) 

Here S stands for the set of quantum numbers that characterize the production 
amplitude A projected on the hyperon-nucleon partial waves. Since we restrict 
ourselves to s-waves in the system, S corresponds directly to the total 
spin of the YN system. In order to simplify the notation we will omit the spin 
index in the following. 

The integrals in Eq. (3) receive contributions from the various possible final- 
state interactions, namely AK^ NK as well as AA^. The first two can be sup- 
pressed by choosing the reaction kinematics properly and thus we may neglect 
them for the moment — we come back to their possible infiuence below — to get, 
for m^ > ml, 
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D{s, t, m^) = A{s, t, m^)e-'^ sin 5, (5) 
where 5 is the elastic AiV (^^o or ^5*1) scattering phase shift ^ . 
The solution of Eq. (3) in the physical region becomes (see Refs. [28,29,30]) 



74(s, t, vrt' 



^u{m^ + iO) 



1 r 



m' ^ — 



(6) 



where, in the absence of bound states. 



u 



TT Jm?. m' ^ — Z 



(7) 



In a large-momentum transfer reaction, the only piece with a strong depen- 
dence on m? is given by the exponential factor in front of the integral in Eq. 
(6). We may thus define 



^4(5, m^) — exp 



1 r"^ 

TT Jml m' 



^ — rn? — iO 



dm! 



1 2 



$(s,i,m'). 



(8) 



where $(s,t,m^) is a slowly varying function of irP'. Henceforth we suppress 
the dependence of the amplitude on s and t in our notation. In the literature 
A{m?) is known as the enhancement factor [25]. 

It is interesting to investigate the form of A{m^) for phase shifts that are given 
by the first two terms in the effective range expansion. 



p' cig{8{m^)) = -1/a + {l/2)rp'^ , (9) 

over the whole energy range. Here p' is the relative momentum of the final 
state particles under consideration in their center of mass system. A{m'^) can 
then be given in closed form as [31] 



^ ^ -l/a+ {r/2)p'^ -ip' ^ ^ ' ^ ^ 



where a = l/r(l + yl — 2r/a). Note that for a ^ r, as is almost realized 
in the ^Sq partial wave of the pp system, the energy dependence of A(m^) is 



^ Obviously this expression holds only for those values of that are below the 
first inelastic threshold. We will ignore this here and come back to the role of the 
inelastic channels later. 
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Fig. 1. The integration contours in the complex m! ^ plane to be used to derive Eqs. 
(12) (a) and (15) (b). The thick lines indicate the branch cut singularities. 

given by 1/(1 + iap') as long as p' <C 1/r. This, however, is identical with the 
energy dependence of pp elastic scattering (if we disregard the effect of the 
Coulomb interaction). Therefore one expects that, at least for small kinetic 
energies, pp elastic scattering and meson production in A^A^ collisions with a 
pp final state exhibit the same energy dependence [25,26], which indeed was 
experimentally confirmed by the measurements of the reaction pp — > ppn^ [27] 
close to threshold. 

The situation is different, however, for interactions where the scattering length 
is of the same order as the effective range, because then the numerator of 
Eq. (10) introduces a non-negligible momentum dependence. This observation 
suggests that a large theoretical uncertainty is to be assigned to the scatter- 
ing length extracted in Ref. [19]. If, furthermore, even higher order terms in 
the effective range expansion are necessary to describe accurately the phase 
shifts, as might be the case, e.g., for the YN interaction [15], no closed form 
expression can be given anymore for A{m'^). Thus one would have to evaluate 
numerically the integral given in Eq. (7) for specific models and then compare 
the result to the data. This procedure is not very transparent. Therefore we 
propose another approach to be outlined in the next subsection. Starting from 
Eq. (8) we demonstrate that the scattering length can be extracted from the 
data directly. Furthermore we show how - at the same time - an estimate of 
the theoretical uncertainty in the extraction can be obtained. 



2.2 Extraction of the scattering length 

Our aim is to establish a formalism that facilitates the extraction of the AN 

scattering length from experimental information on the invariant mass dis- 
tribution. To do so we derive a dispersion relation for |A(m^)p. We will use 
a method similar to the one utilized in Refs. [32,33]. First, notice that - by 
construction - the function 
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log {Ainfy^inf)} ^ 1 1 r '^(^'^^ dm''' (11) 
- ml yjm^ - ml^ '< m'^-m'^-iO 

has no singularities on the physical sheet except for the cut from m^ to infinity. 
In addition, its value below the cut is equal to the negative of the complex 
conjugate of the corresponding value above the cut. Therefore, from Cauchy's 
theorem we get 



/ 2^ 



log{.4(//r)/<I'(///2)} 1 log|.4(///'2)/<I.(„/2)i2 



r^2m^ini3!lJL_dm!\ (12) 



^Jm^ — ml 27ri Jmg ^777,' 2 _ ^2 j-^/ 2 _ ^2 _ ^q^ 

Taking the imaginary parts of Eqs. (11) and (12) one gets 



5{m^) -D log |A(m'^)/$(m' 



/2M2 



----Pi -"^i"^-- ^'^^-^ ^1 dm'^ . (13) 

Jm^ — ml ^"^ "^"*o ym'^ — mo (m'2 — 777,2) 

Note that for small AA^ kinetic energies we may write m' = ml + ''^oP ^ 

where ji — j^^^^^^. Thus, for small invariant masses, ^Jm^ — ml oc p', 
which clearly demonstrates that the left hand side of Eq. (13) converges to 
the (negative of the) scattering length times ji/mo for m^ — > ml- 



The integral in Eq. (13) depends on the behaviour of A{m''') for large m'^ 
which is not accessible experimentally. Therefore, in practice the integration 
in Eq. (13) has to be restricted to a finite range. In fact, a truncation of the 
integrals is required anyway because in our formahsm we do not take into ac- 
count exphcitly (and we do not want to) the complexities arising from further 
right-hand cuts caused by the opening of the SA^, nAN, etc., channels. Thus, 
let us go back to the definition of A in Eq. (8). A significant contribution of 
the integral there stems from large values of m' ^. Those contributions depend 
only weakly on m^, at least for values in the region close to threshold. 
Therefore, they can be also absorbed into the function $ defined in Eq. (8). 
Consequently, we can write 



A(m^) = exp 



X f^max 



S{m' 



1 2^ 



TT Jml m 



I 2 



777.2 _ 



-dm 



I 2 



(14) 



where ^(m^^^, 777^) = ^{m^)^^i^^^{m') with 



""max 



1^777^) = exp 



-/ 



5(777' 2) 



2,„, m''^-m^ - 70' 



-dnC 



1 2 



(15) 
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The quantity ni^yiax be chosen by physical arguments in such a way 

that both $(m^) and $m2^^^(m^) vary slowly over the interval (^^^O' ""^maa;)- 
particular, it is crucial that it can be chosen to be smaller than the invariant 
mass corresponding to the TiN threshold. Note also that, strictly speaking, 
Eq. (15) is not valid anymore in the presence of inelastic channels. Rather one 
should view Eq. (14) as the actual definition of But we will still use Eq. 
(15) later to estimate the uncertainty of Eq. (2). 

We should mention that the integral in Eq. (14) contains an unphysical sin- 
gularity of the type log (m^^^. — w?) - which is canceled by a corresponding 
singularity in $(m^„3.,m^) - but this again does not affect the region near 
threshold. 

The method of reconstruction of 5{w?) is similar to the one used for the infinite 
range integration. The function 



log{^(m')/l>(m^,„m=^)} 
^(m2 -mg)(m^„^ -m2) 



ml){ml 



X 



1 

TT Jm 



I 2 



^0 



12 



(16) 



2 

max • 



has no singularities in the complex plane except the cut from rn^ to m. 
Again, its value below the cut is equal to the negative of the complex conjugate 
of its value above the cut. Hence, applying Cauchy's theorem, one gets 



log{A(m^)/l>(m^, ^,m=^)} 
— m?) 



{w? — mQ){m: 

1 
2wi 



2 

max 



log|A(m'2)/$(mL^,m'2)|2 



{m' 2 — mo) {iTi'^ax ~ ^) ("^' ^ ~ ~ ^0) 



dm' 



I 2 



and, accordingly. 



\jm^ — rr 

--pf 

2tT Jn, 



^2 

"max 



/2\\2 



\og\A(m'^)/^ml^^,m 
\Jm''^ — mg {m' ^ — w?) \ 



mt, 



^max-^'^ 



dm 



I 2 



(17) 



It is important to stress that the principal value integral vanishes for a con- 
stant argument of the logarithm as long as m^ < m^ < m"^^^. Therefore, if 
the function $ depends only weakly on m' ^ — as it should in large-momentum 
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transfer reactions — it can actually be omitted in the above equation. In addi- 
tion, 



dm' "^dt 



where as denotes the partial cross section corresponding to As - and where we 
included again the explicit S dependence as a reminder. Thus we can replace 
|A(m^)P in Eq. (17) by the cross section, where all constant prefactors can 
be omitted again because the integral vanishes for a constant argument of the 
logarithm. The result after performing these manipulations is given in Eq. (2). 

Since Eq. (17) gives an integral expression not only for the threshold value 
of the phase shift but also for the energy dependence, one might ask whether 
it is possible to extract, besides the scattering length, also the effective range 
directly from production data. Unfortunately, this is not very practical. By 
taking the derivative with respect to of both sides of Eq. (17) one sees that 
one can get only an integral representation for the product a|((2/3)a5 — r^) 
but not for the effective range rs alone. Thus, although the corresponding 
integral is even better behaved than the one for as in Eq. (2) (it has a weaker 
dependence on m^^^), the attainable accuracy of rs will always be limited by 
twice the relative error on a^. In fact, a more promising strategy to pin down 
the YN low energy parameters might be to fix the scattering lengths from 
production reactions and then use the existing data on elastic scattering to 
determine the effective range parameters. In this case one would benefit from 
having to do an interpolation between the results at threshold, constrained 
from production reactions, and the elastic YN data available only at somewhat 
higher energies. This should significantly improve the accuracy as compared 
to analyses as, e.g., in Ref. [8] which rely on elastic scattering data alone. 



2.3 Relevant scales and error estimates 



The next step is to estimate an appropriate range of values for m^^^. as well 
as the systematic uncertainty of the method. Naturally, m^^^ has to be large 
enough to allow resolution of the relevant structures of the final state inter- 
action. (Note that rn^ leads to a vanishing value for the scattering 
length, c.f. Eq. (2).) On the other hand it should be as small as possible in order 
to ensure that the h.N system is still produced predominantly in an jS-wave 
but also, as mentioned, to avoid the explicit inclusion of further right-hand 
cuts such as the one resulting from the opening of the TiN channel. 

Eq. (10) suggests that we should choose m^^^, large enough to allow values of 
p' such that ap' ~ 1. Thus, for non-relativistic kinematics, we find m^^^. ~ 
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Fig. 2. Dependence of the extracted scattering lengths on the value of the upper 
limit of integration, emax- Shown is the difference from the exact results for the 
spin singlet scattering length (left panel) and the spin triplet scattering length at 
(right panel). The solid and the dot-dashed line correspond respectively to model 
NSC97a and NSC97f of Ref. [5] and the dashed one corresponds to the YN model 
of Ref. [6]. The shaded area indicates the estimated error of the proposed method 
and the arrows indicate the value for e^nax as estimated based on scale arguments. 



mo(l + l/(/imoa^)). A more transparent quantity is emax defined by 



rrimax - mo 



i.e. the maximum kinetic energy of the AA'^ system for which it contributes 
to the integral of Eq. (2). For a ~ 1 fm we get e^ax ~ 40 MeV. This is well 

below the UN threshold which corresponds to an energy of 75 MeV and the 
threshold for the ttAA^ channel which is at 140 MeV. 



Let us now estimate the uncertainty. Except for the neglect of the kaon-baryon 
interactions, Eq. (17) is exact. Therefore 5a'^^^'> - the theoretical uncertainty 
of the scattering length extracted using Eq. (2) - is given by the integral 
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Sa^^- lim J_(^^ + ^A 



f 

Jm 



drn'^ (18) 



Since log |<l>(m^„^, m^)^ = log|<l>(m2)|2 + log |<l>„2^^^(m2)|2, we may write 
5a*^*^^ = 5a^^^'^^ + ^a"^""*^, where the former, determined by $(m^), is con- 
trolled by the left hand cuts and the latter, determined by $TO2^^^(m^), by the 
large energy behavior of the KN scattering phase shifts. The closest left hand 
singularity is that introduced by the exchange of light mesons (tt and K) in 
the production operator and it is governed by the momentum transfer. Up 
to an irrelevant overall constant, we may therefore estimate the variation of 
$ ~ 1 + /s;(p7p)^) where we assume k to be of the order of 1. Evaluation of 
the integral (18) then gives 

^aC'^^) ^ ^^{p'maJp') - 0.05 fm . 

Here we used p'^ax ^ = '^l^^max with Cmax ~ 40 MeV and p ~ 900 MeV, where 
the latter is the center-of-mass momentum in the initial NN or 'yd state that 
corresponds to the KAN threshold energy. Concerning ^a™™"^ we start from 
the definition of $^2 (m^) in Eq. (15) from which one easily derives 



\Sa" 



6{y)dy 



2V3/2) 



< 



(19) 



(m 



mz 



2 

max 



itlq). Thus, in order to obtain an estimate 



where y"^ 

for 50™"™""= , wc need to make an assumption about the maximum value of the 
elastic AA^ phase for > m^^^. In addition, it is important to note that 
the denominator in the integral appearing in Eq. (19) strongly suppresses 
large values of y. Since for none of the considered YN models [2,3,4,5,6] does 
Smax exceed 0.4 rad, we arrive at the estimate 5a'"'"''"' ~ 0.2 fm. Note that 
when using the phase shifts as given by those models directly in the integral, 
the value for 50™™""= is significantly smaller, since for all models considered the 
phase changes sign at energies above m'^^^. Combining the two error estimates, 
we conclude 

Sa < 0.3 fm . 



Another source of uncertainty comes from possible final state interactions in 
the KN and KA systems. Their effects have been neglected in our consid- 
erations so far for a good reason: their influence is to be energy dependent 
and therefore no general estimate for the error induced by them can be given. 
For example, baryon resonances are expected to have a significant infiuence 
on the production cross sections [36]. However, it is possible to examine the 
influence of the meson-baryon interactions experimentally, namely through 
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a Dalitz plot analysis. This allows to see directly, if the area where the AN 
interaction is dominant is isolated or is overlapping with resonance structures 
[37] . In addition, to quantify the possible influence of the meson-baryon inter- 
actions on the resulting scattering lengths, the extraction procedure proposed 
above is to be performed at two or more different beam energies. If the ob- 
tained scattering lengths agree, we can take this as a verification that the 
result is not distorted by the interactions in the other subsystems because 
their influence necessarily depends on the total energy, whereas the extraction 
procedure does not. 



3 Test of the method: comparison to model calculations 

The most obvious way to test Eq. (2) would be to apply it to reactions where 
the involved scattering lengths are known as it is the case for any large mo- 
mentum transfer reaction with a two nucleon final state. Unfortunately, as far 
as we know there are is no experimental information on nn or pn mass spec- 
tra with sufficient resolution to allow the application of Eq. (2). For existing 
invariant mass spectra with a pp final state, as given, e.g., in Ref. [38] for 
the reaction pp ppn^, our formula is not applicable due to the presence of 
the Coulomb interaction that strongly distorts the invariant mass spectrum 
exactly in the regime of interest. Moreover, one has to keep in mind that the 
authors of Ref. [38] were forced to assume already a particular dependence 
of the outgoing pp system at small values of for their analysis because of 
the limited detector acceptance. Thus, a test with those data would not be 
conclusive, anyway. 

Therefore, we chose a different strategy in order to test our method. We apply 
Eq. (2) to the results of a microscopic model for the reaction pp K~^pA 
[35]. (We expect that application to the reaction ^fd K^nA gives similar 
results.) In this model the K production is described by the so-called tt and K 
exchange mechanisms. Thereby a meson (tt or K) is first produced from one 
of the nucleons and then rescattered on the other nucleon before it is emitted 
(cf. Ref. [35] for details). The rescattering amplitudes {-kN KY, KN 
KN) are parameterized by their on-shell values at threshold. The YN final 
state interaction, however, is included without approximation. Specifically, 
the coupling of the AA^ to the UN system is treated with its full complexity. 
For the present test calculations we used all the YN models of Ref. [5] as 
well as the YN model of Ref. [6]. Thus, the test calculation will allow us to 
see a possible infiuence of the EA^ channel as well as that of the left hand 
singularities induced by the production operator. In addition, we can study 
the effect of the value of the upper integration limit m^^^ or, equivalently, 
^max on the scattering length extracted. 
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2060 2080 2100 



nip^ [MeV] 

Fig. 3. pp K^X inclusive missing mass spectrum at Tp = 2.3 Gev for kaon labo- 
ratory scattering ang les = 10±2° [39]. The solid line is our fit as described in the 
text. Note that the scattering length extracted is insensitive to the parameterization 
of the data used (see Appendix for details). 

We now proceed as follows: we calculate the invariant mass spectra needed 
as input for Eq. (2), utilizing the YN models mentioned above. Then we use 
that equation to extract the scattering lengths from the invariant mass plots. 
In Fig. 2 the difference between the extracted scattering length and the exact 
one is shown as a function of emax for several YN models. The direction of 
kaon emission was chosen as 90 degees in the center of mass. Considered are 
the two "extreme" models of Ref. [5], namely NSC97a, with singlet scattering 
length tts = —0.71 fm and triplet scattering length at = —2.18 fm, and NSC97f 
[as = —2.51 fm and at = —1.75 fm), and the new Jiilich model (a^ = —1.02 
fm and at = —1.89 fm). The last was included here because it should have 
a significantly different short range behavior compared to the other two. As 
can be seen from the figure, for emax — 40 MeV in all cases the extracted 
value for the scattering length does not deviate from the exact one by more 
than the previously estimated 0.3 fm, as is indicted by the shaded area in 
Fig. 2. This is in accordance with the error estimates given in the previous 
section, where emax ~ 40 MeV - in the figure highlighted by the arrows - 
was deduced from quite general arguments. In addition, this investigation also 
suggests that the T,N cuts, included in the model without approximation, do 
not play a significant role in the determination of the scattering length. 
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4 Application to existing data 



As was stressed in the derivation, Eq. (2) can be applied only to observables 
(cross sections) that are dominated by a single partial wave. This is indeed 
the case for specific polarization observables for near-threshold kinematics 
as we show in the appendix B for the reaction pp NKA. Unfortunately, 
at present no suitable polarization data are available that could be used for 
a practical application of our formalism. Thus, in order to demonstrate the 
potential of the proposed method, we will apply it to the unpolarized data 
for the reaction pp — > K'^X measured at the SPES4 facility in Saclay [39]. 
Obviously, since we only consider the small invariant mass tail of the double 
differential cross section, we automatically project on the pA S-waves. The 
data, however, still contains contributions from both possible final spin states 
- spin triplet as well as spin singlet. Thus the scattering length extracted 
from those data is an effective value that averages over the spin-dependence 
of the pA interaction in the final state and also over the spin dependence of 
the production mechanism with unknown weightings. It does not allow any 
concrete conclusions on the elementary pA (}So or ^5*1) scattering lengths. We 
should mention, however, that some model calculations of the reaction pp 
K^pA [35,40] indicate that the production mechanism could be dominated by 
the spin-triplet configuration. If this is indeed the case then the production 
process would act like a spin filter and the scattering length extracted from 
the Saclay data could then be close to the one for the pA '^'Si partial wave. 

For our exemplary calculation we use specifically the data on the pp K^X 
inclusive missing mass spectrum at Tp = 2.3 GeV for kaon laboratory scat- 
tering angles = 10±2°^ [39] (cf. Fig. 3). For convenience we represent 
those data in terms of simple analytical functions (the result of the best fit is 
shown as the solid line in the figure) as described in appendix A. Please note, 
however, that the advantage of the proposed method is that the scattering 
length extracted is independent of the particular analytical parameterization 
employed; one can even use the data directly. In particular we do not need 
to assume that the elastic AA^ interaction can be represented by the first two 
terms of an effective range expansion. The analytical representation of the 
data is then used to evaluate the corresponding scattering length utilizing Eq. 
(2). Details of the fitting as well as how we extract the uncertainty in the 
scattering length can be found in the appendix. 

We find a scattering length of (—1.5 ±0.15 ±0.3) fm, where the first er- 
ror corresponds to the uncertainty from the data and the second one is the 
theoretical uncertainty discussed previously. Thus, already in the case of an 



^ For the given energy and ernax='^^ MeV a fixed kaon angle is almost equivalent 
to a fixed t as is required for the dispersion integral. 
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experimental resolution of about 2 MeV the experimental error is smaller than 
the theoretical one. 



5 Summciry 

In this paper we presented a formalism that allows one to relate spectra from 
large-momentum transfer reactions, such as pp — > K~^pA or 70? — > K'^nA, 
directly to the scattering length of the interaction of the final state particles. 
We estimated the theoretical error of the analysis to be less than 0.3 fm. 
This estimate was confirmed by comparing results obtained with the proposed 
formalism to those of microscopic model calculations for the specific reaction 
pp K^Ap. 

The formalism can be applied only to observables (cross sections) that are 
dominated by a single partial wave. This requirement is fulfilled by specific 
polarization observables for AN invariant masses near threshold as we show 
in the appendix B for the reaction pp — > NKA. Corresponding experiments 
require a polarized beam or target in order to pin down the spin-triplet scatter- 
ing length and polarized beam and target for determining also the spin-singlet 
scattering length. 

Since at present no suitable polarization data are available we demonstrated 
the potential of the proposed method by applying it to unpolarized data 
from the reaction pp — >■ K^X measured at the SPES4 facility in Saclay [39]. 
Thereby a scattering length of (—1.5 ± 0.15 ± 0.3) fm was extracted. This 
more academic exercise demonstrates that the resolution of 2 MeV of the 
SPES4 spectrometer is already sufficient to obtain scattering lengths with an 
experimental error of less than 0.2 fm. Such an accuracy would be already an 
improvement over the present situation. 

The method of extraction of the YN scattering lengths discussed here should 
be viewed as an alternative method to proposed analyses of other reaction 
channels involving the YN system like K~d — > ^nA. The assumptions that go 
into the analyses are very different and therefore carrying out both analyses 
is useful in order to explore the systematic uncertainties. 
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A Method of the Analysis 



In order to evaluate the integral in Eq. (2) for the experimental results of 
Ref. [39] we first fit the data utilizing the following parameterization for the 
amplitude squared (c.f. Eq. (17)): 



\A{m) 



exp 



(A.l) 



The advantage of such a parameterization is that the integral in Eq. (2) can 
be calculated analytically and one gets for the scattering length 



1 



a{C,,C^)^-^-Cl. 

2 \ ytriisffnA^ 



rrin 



[mt 



c|)K-c|)3- 



(A.2) 



One can extend Eq. (A.l) easily by including in the exponent additional terms 
such as C|/ (m^ — C|) if required for a satisfactory analytical representation of 
the data. The accuracy of the presently available experimental results, how- 
ever, did not necessitate such an extension, for already the form given in Eq. 
(A.l) ensured a per degree of freedom of less than 1. 

In Ref. [41] it is stressed that there is a calibration uncertainty in the data 
of Ref. [39]. The actual value of this uncertainty was determined by including 
as a free parameter a shift Am in the fitting procedure. In this way a value 
of Am = 1.17 MeV was found. We use the same value. In addition, the finite 
resolution of the detector cr^ = 2 MeV has to be taken into account. Thus the 
function 



rf^(j(m) 
dmdflK 



'(fa{m') 



exp 



— {m — m' — Amy 



dm' ,(A.3) 



is to be compared to the data, where {d'^a{m')/{dm'dVtK))o is calculated from 
the amplitude Eq. (A.l). 

The probability for the data to occur under the assumption that the cross 
section given in Eq. (A. 3) is true is given by the Likelihood function. Here we 
assume the data distribution to be Gaussian: 



i:({data}|C70, C71, C2; /) oc exp 



-\x\Co,C,,C2) 



(A.4) 



where the standard function appears in the exponent. The letter I that 
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appears in the argument of L is meant to remind the reader that some as- 
sumptions had to be made in order to write down Eq. (A.4). 

The distribution of interest to us is the probabihty distribution of the scatter- 
ing length £(a|{data}; /), given the data, that can be written as [42] 



£(a|{data};/) = 

H f lldCiC{{dsitsi}\C0,Cl,C2; I)6{a - a{Ci,C2)) ,(A.5) 

•' i 

where is a normahzation constant to be fixed through / C{a)da — 1 and 
a(Ci,C2) is given in Eq. (A.2). 

It turns out that the LikeUhood function of Eq. (A.4) is strongly peaked. 
For the numerical evaluation of Eq. (A. 5) we therefore linearize the individual 
terms with respect to the parameters Cj. Within this linear approximation the 
resulting probability density function C{a) also takes a Gaussian form. To test 
that there is no significant dependence of the result on the parameterization 
used for the data we also used the Jost function in the analysis. The results 
for the scattering lengths determined with the two different methods agreed 
within the statistical uncertainty. 



B Spin dependence of the reaction pp — > ApK'^ 



The aim of this appendix is to show how polarization measurements for the 
reaction pp — > ApK'^ can be used to disentangle the spin dependence of the 
production cross section. 

In terms of the so called Cartesian polarization observables, the spin-dependent 
cross section can be written as [43] 



+ J2{PbUPt)jM0 + ... 



(B.l) 



where cro(^) is the unpolarized differential cross section, the labels i,j and k 

can be either x, y or z, and P^, Pt and Pf denote the polarization vector of 
beam, target and one of the final state particles, respectively. All quantities are 
functions of the 5 dimensional phase space, abbreviated as ^. The observables 
shown explicitly in Eq. (B.l) include the beam analysing powers A^q, the 
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corresponding quantities for the final state polarization DQi, and the spin 
correlation coefficients A^j . Polarization observables not relevant for this work 
are not shown explicitly. All those observables that can be defined by just 
exchanging and Pt, such as the target analyzing power Aoi, are not shown 
explicitly. 

The most general form of the transition matrix element may be written as [37] 



M = H{II') + iQ ■ {SI') + lA ■ {S'X) + {Si Sj ')Bij , (B.2) 

where S — (x2'^y^Xi) ^^'^ ^' — (x4'^'^2/(xD^) ^e used for spin triplet 

initial and final states, respectively, and I — (xl^^yXi) ^^'i ^' — (x4'^j/(X3 )^) 
are to be used for the corresponding spin singlet states. Here it is assumed 
that the outgoing baryons are non relativistic, as is the case for the kinematics 
considered. 

Given this form it is straight forward to evaluate the expression for the various 
observables. E.g. 
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4(70 = + \Qm? + \Am? + \B,nu? . (B.3) 

AAoiGQ = ieikn {QlQn + Bl^Bnrn) + 2Im - Q*//) , (B.4) 

ADoiao = -ieikn {AIA„ + B*^,B^^) - 2Im {B*,,Q„^ - A*H) , (B.5) 

4Ai<70 - -W + \Qm\' - \Am\' + \Bmn\' ' 2\Qi\' - 2| Ami' , (B.6) 

where the indices m, k and n are to be summed. By definition, the spin 
triplet final state amplitudes contribute to A and B only, whereas the spin 
singlet amplitudes contribute to H and Q. By definition, the spin triplet states 
interfere with the spin singlet states only, if the final state polarization is 
measured. In the examples given above this is the case for D^i only. This 
simplifies the analysis considerably. Eq. (2) needs as input only those cross 
sections that have low relative energies, where we can safely assume the AN 
system to be in an S'-wave. Thus, in order to construct the most general 
transition amplitude under this condition, we need to express the quantities 

-* — * — * 

H, Q, A as well as B^j in terms of the vectors k (the momentum of the 
kaon) and p (the momentum of the initial protons) — c.f. Fig. B.l. The Pauli 

— * 

Principle demands for a proton-proton initial state, that S appears with odd 

powers of p and terms that do not contain S appear with even powers of p. 
Correspondingly, due to parity conservation these amplitudes need to be even 
and odd in k respectively, since the kaon has a negative intrinsic parity. We 
may therefore write 

H = 0, 

Q = siik ■ p)k + 
A = tip{p-k) + t2k, 
Bij = Sijk (t^kkik ■ p) + Upi^ + {k X p)i {t^kj +h{p- k)pjj , (B.7) 

where p — p/\p\- The amplitudes Sj and ti appearing in the above equations 
are functions of p^, /c^ and {k ■ p)"^. Now we want to identify observables that 
either depend only on (some of) the Sj and thus are sensitive to the spin singlet 
AN interaction only or on (some of) the ti and thus are sensitive to the spin 
triplet AA^ interaction only. It should be stressed in this context that it is 
entirely due to the relation H — - a, direct consequence of the selection rules 
- that one can disentangle the spin triplet and the spin singlet channel without 
making assumptions about the kaon partial wave, as was first observed in Ref. 
[44]. 

Using Eqs. (B.7) all observables can be expressed in terms of the amplitudes 
Si and ti, where the former correspond to a spin singlet YN pair in the final 
state, while the latter correspond to a spin triplet pair. We find 

AoiCTo ^^{kx P)i (a + I3{k ■ p) + -i{k ■ pf) (B.8) 
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where 



a = 1111(^:^2 + ^5^2^; , (B.9) 
f3^1m{tlh + s*si) , (B.IO) 
-f^lm{tlh + tlh + tlt2-tlh) . (B.ll) 

We may express the results in the coordinate system defined by the beam 
along the z axis and the x and y directions through the polarizations (c.f. Fig. 
B.l). Then we find for the analysing power 



^oyO'o = —-k'^P sin(26') cos(0) 

+ sin(^) cos(0) (spin triplet only) , (B.12) 

— * — * 

where we used (p x k)y = /c sin(6') cos(0) and p ■ k — kcos{9). Thus, the 
contributions from all those partial waves where the final AA^ system is in 
an ^5*0 state vanish when the kaon goes out in the xy-plain. In addition, 
the term proportional to j3 is the only one that is odd in cos(^^), and thus 
any integration with respect to the angle 9 of AoyaQ symmetric around tt/2 
removes any spin singlet contribution from the observable. Please note that for 
any given energy it should be checked whether the range of angular integration 
performed is consistent with the requirement that the reaction was dominated 
by a large momentum transfer - kaon emission at 90 degrees maximizes the 
momentum transfer and thus minimizes the error in the extraction method. 

Analogously one can show for the combination (1 — Axx)cro that 

(1 - Axx)ao = i + ig^l' + \Bxm\^) 
I 

= -k^\si\^sm^{29) cos^(0) + (spin triplet only) , (B.13) 

8 

i.e. that it also allows to isolate the amplitudes with the YN system in the 
spin triplet. Furthermore, 



(1 + + Ayy — Azz)(^Q 

— \Q z\ ^" \Bzm\ 

= ^k^\tz\^ sin^(2^) + (spin singlet only) (B.14) 

and therefore allows to isolate the spin singlet amplitudes. (Note that in both 
cases a summation over m is to be performed.) Thus, for all those partial 
waves where the final YN system is in the ^Sq final state, (1 — A^x) vanishes. 
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when the kaon goes out in the yz- or in the xy-plain and for all those partial 
waves where the final YN system is in the ^Si final state, {1 + A^.^.^ Ayy — A^z) 
vanishes, when the kaon goes out in the x?/-plain or in forward direction. These 
results hold independent of the orbital angular momentum of the kaon. 

So far we assumed that higher partial waves of the AiV system do not play 
a role - since we restricted ourselves to small relative energies. However, it is 
also possible to check this experimentally. First of all the angular distribution 
of the YN system da/dQp/, where p' is the relative momentum of the two 
outgoing baryons, needs to be fiat. However, in the presence of spins even fiat 
angular distributions can stem from higher partial waves [45]. To exclude this 
possibility the A polarization can be used. For this test we need an observable 
that vanishes in the absence of higher AA^ partial waves. The easiest choice 
here is Dqi. It is straight forward to show that 

Doi (x{px k)i , 

as long as there are only S waves in the AA^ system. Therefore, in the absence 
of higher partial waves in the YN system Dq^ has to vanish. 
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